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, ABSTRACT 

■ We present near infrared J-, H- and K-band photometry and optical spectroscopy 

of low-mass star and brown dwarf (BD) candidates in the Pleiades and Praesepe 

{SJ . open clusters. We flag non-members from their position in K, I-K and J, J-K colour- 

J> ' magnitude diagrams (CMDs), and J-H, H-K two-colour diagrams. In general, the 

theoretical NextGen isochrones of the Lyon Group fit the K, I-K CMDs well for stars 

■ with I-K^l.5-3.5. However Pleiades stars with K~10.5-13 (Mk —5-7.5) are rather 

redder than the isochrones. We also identify this effect amongst aPer sources from the 
literature, but find no evidence of it for field stars from the literature. The NextGen 
isochrones fit the J, J-K CMDs of both clusters very well in this photometric range. It 
is possible that the I-K colour of youthful stars is affected by the presence of magnetic 
activity. The Lyon Group's Dusty isochrones fit both K, I-K and K, J-K Pleiades 

■ CMDs well for I-K~4.3-6/J-K~l. 1-1.4. In between these colour ranges the Pleiades 

cluster sequence comprises three portions. Starting at the bluer side, there is a gap 

q ■ where very few sources are found (the gap size is AI~0.5, AJ^AK^0.3), probably 

$H ' resulting from a sharp local drop in the magnitude-mass relation. Then the sequence 

, is quite flat from I-K~3.5-4. Finally, the sequence turns over and drops down to join 

■ the Dusty isochrone. We also compare model atmosphere colours to the two-colour 

diagrams of the clusters. The NextGen models are seen to be ~0.1 too blue in H-K, 

1— 1 , and ~0.1 too red in J-H for T e ff >4000K. However, they are in reasonable agreement 

>^ ■ with the data at T e// -3200K. For T e// ~2800-3150K, the colour of Pleiades and 

Praesepe sources are significantly different, where Praesepe sources are ~0.1 bluer in 
J-H and up to ~0.1 redder in H-K. These differences could result from gravity sensitive 
molecular opacities. Cooler Praesepe sources then agree well with the dusty models, 
suggesting that dust is beginning to form in Praesepe sources around 2500K. However, 
Pleiades sources remain consistent with the NextGen models (and inconsistent with 
the dusty models) down to T e //s of ^2000K. It is possible that dust formation does 
not begin until lower T e y/s in sources with lower surface gravities (and hence lower 
atmospheric pressures). We also identify unresolved binaries in both clusters, and 
estimate mass ratios (q) for Pleiades BDs. Most of these have q>0.7, however, 3/18 
appear to have lower q values. We determine the binary fraction (BF) for numerous 
mass ranges in each cluster, and find that it is generally rising towards lower masses. 
We find a BD BF of 50^q%. We also find some evidence suggesting that the BF-q 
distribution is flat for 0.5-0.35M Q , in contrast to solar type stars. 

Key words: stars: low-mass brown dwarfs - open clusters and associations: individ- 
ual: Pleiades Praesepe 
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1 INTRODUCTION 

Low-mass star and brown dwarf (BD) populations in open 
clusters make excellent test-beds for our theoretical under- 
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standing of these objects. With a well constrained cluster 
age, and assuming uniform composition, we can effectively 
compare theoretical isochrones to observed colour magni- 
tude diagrams (CMDs) and two-colour diagrams without 
the complications introduced by a spread in metallicity and 
a range in surface gravities for a given mass. Furthermore, 
for most open clusters, the depth effect is small, and the 
single star sequence is therefore quite tight on CMDs. ft is 
thus possible to study unresolved binarity amongst cluster 
members. 

The T e f f range of young to intermediate age low- mass 
stars and BDs corresponds to M and L spectral types, and 
modeling the atmospheres of such objects is very compli- 
cated. The convection zone penetrates deeply into the opti- 
cally thin outer layers due to H2 dissociation (Burrows et al. 
1993). Molecular opacities such as TiO and H2O comprise 
millions of lines that must be accurately modeled (Allard 
et al. 1997). Also, dust begins to form in atmospheres at 
around 2800K (Allard et al. 2001), which depletes opacities 
as well as affecting the atmospheric temperature gradient. 
There are many different varieties of dust (Sharp & Huebner 
1990), and their formation rates and grain size distributions 
will depend not only on temperature, but on surface gravity 
and metallicity. Despite these complexities, much progress 
has been made in modeling these atmospheres. The Lyon 
Group has combined structural models (using interior struc- 
ture physics) with model atmospheres (Hauschildt, Allard 
& Baron 1999; Allard et al. 2001) to calculate evolutionary 
models for low-mass stars and BDs, both with and without 
atmospheric dust condensates (Baraffe et al. 1998, hereafter 
NextGen; Chabrier et al. 2000, hereafter Dusty). Calcula- 
tions by other groups include those of Burrows et al. (1997), 
and D'Antona & Mazzitelli (1994). However, only the Lyon 
group models are based on consistent non-grey evolutionary 
calculations relating the mass and age of an object to its 
observational properties (colour and magnitude). 

The Pleiades and Praesepe are two nearby, rich open 
clusters with near solar composition (Crawford & Perry 
1976; Hambly et al. 1995b). With distances of ~125pc and 
160pc respectively, both clusters have small depth effects, 
and are fairly compact on the sky (a few degrees across). 
The Pleiades is ~125Myrs old (Stauffer, Schultz & Kirk- 
patrick 1998), and Praesepe is ~0.6-1.4Gyrs old (Hambly 
et al. 1995b). As a consequence, Pleiades low- mass stars 
(ZeO.3M ) and BDs are not fully contracted, and will have 
lower surface gravities than their equal T e / / counterparts in 
Praesepe. 

In this paper, we derive membership criteria of Pleiades 
and Praesepe low-mass star and BD candidates using near 
infrared photometry and optical spectroscopy. We compare 
NextGen and Dusty theoretical isochrones to our observed 
CMDs and two-colour diagrams, and identify photometric 
ranges where the isochrones agree well with observations. We 
then join these isochrone sections together using additional 
empirical points, and discuss the overall form of the cluster 
sequences. We also use our CMDs and two-colour diagrams 
to identify unresolved binaries. We derive binary fractions 
over several mass ranges in each cluster, and estimate mass 
ratios for the Pleiades BD binaries. Finally, we discuss future 
work. 



2 OBSERVATIONS & DATA REDUCTION 

2.1 The Pleiades & Praesepe samples 

Our Pleiades and Praesepe samples were taken from several 
sources. For the Pleiades, the highest mass sources ((I<13.5) 
were taken from Hambly, Hawkins & Jameson (1993; HHJ 
hereafter). Fainter sources were taken from, the ITP survey 
(Zapatero-Osorio et al. 1999), the CFHT survey (Bouvier et 
al. 1998), the Burrell Schmidt survey (Pinfield et al. 2000; 
BPL hereafter) and the INT WFC survey (Dobbie et al. 
2002b). The fainter surveys cover a combined area of 7.6 
square degrees, and range in photometric depth from 1=19.5- 
21.8 (unless otherwise stated, I refers to Cousins I). These 
fainter surveys are summarised in more detail in Jameson et 
al. (2002). 

Our Praesepe sample was taken from four surveys. 
Proper motion members from the 19 square degree sur- 
vey of Hambly et al. (1995a, hereafter HSHJ) represent our 
highest mass sources (1=10-17.5). Lower mass sources come 
from Pinfield et al. (1997), Pinfield (1997) and Magazzu et 
al. (1998). The Pinfield surveys covered 1 and 6 square de- 
grees, and reached 1=21 and 19.5 respectively. The Magazzu 
survey covered 800 square arcminutes down to 1=21.2, and 
identified 1 BD candidate which they called Roque Pr 1. 

2.2 Photometry 

For the brighter BPL and HSHJ candidates we obtained J- 
H- and K-band magnitudes from the 2MASS 2nd incremen- 
tal point-source catalogue at the NASA/IPAC InfraRed Sci- 
ence Archive (Skrutskie et al. 1995). We used a search radius 
of 3 arcseconds to allow for positional uncertainties. For the 
ITP, CFHT and INT WFC surveys, we obtained photometry 
from the literature for as many sources as was available (see 
Jameson et al. 2002). Most of the fainter BPL candidates 
had K-band measurements from Pinfield et al. (2000), and 
several Praesepe candidates had K-band measurements from 
Hodgkin et al. (1999). For the remainder of these samples 
we measured J-, H- and K-band photometry using UFTI on 
the United Kingdom Infrared Telescope (UKIRT) on Mauna 
Kea, Hawaii, from 3-5 January 2000, and in service time on 
8 October 2001. The sky conditions were photometric. Our 
observing strategy for candidates without a K-band mag- 
nitude was to observe them first at K, and only go on to 
measure J- and H-band magnitudes if the targets had a suf- 
ficiently red I-K colour to be consistent with membership of 
the cluster. For the remaining candidates, we simply mea- 
sured J- and H-band magnitudes, as well as K-band if the 
previous K uncertainty was >0.1. Each observation com- 
prised a five-point dither pattern, with a windowed 512 x 
512 pixel readout. Total integration times at each of J, H 
and K were 150s, 300s and 540s, for objects with 1=17- 
18.5, 1=18.5-20 and I>20 respectively. The raw images were 
de-biased, dark subtracted, flat-fielded, and combined into 
mosaics using ORACDR. Aperture photometry was then 
extracted using GAIA. As the UKIRT photometry was ob- 
tained with the "Mauna Kea Observatory" (MKO) filter set 
(Simons & Tokunaga 2002), the photometry from the lit- 
erature was transformed into the MKO system using the 
transforms of Hawarden et al. (2001) and Carpenter (2001) 
(old UKIRT system and 2MASS respectively). The photom- 
etry is given in the appendices. Where there were multiple 
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Table 1. Spectral indices and equivalent widths of the Pleiades BD candidates. 



Source 


SpT 


PC3 


A 


VO 


13 


12 


H Q EW 


KI 7655 


KI 7699 


Nal 8183+8195 
















(A) 


EW (A) 


EW (A) 


EW (A) 


BPL 327 


dM7.1 


1.55 


1.52 


2.45 


2.33 


1.03 


-5±2 


-12±3 


-10±3 


-7±1 


BPL 45 


dM8.1 


1.76 


1.69 














-7±2 



measurements available of an object, we selected those with 
small errors (generally trying to avoid measurements with 
errors larger than ~0.1), and took an average. 

Our optical photometry was taken from the literature. 
All I-band photometry was transformed into the Cousins 
system. Photographic In magnitudes from HSHJ were trans- 
formed using Bessell (1986). To transform the Ikp magni- 
tudes of Pinfield (1997) and Pinfield et al. (2000) we de- 
rived the following relation; I c -Ik p=0.049(Ik p-K)+0.069 
(for 2.3 < Ic — K < 5.8). This relation was obtained us- 
ing synthetic photometry derived from spectra, using the 
method and data from Dobbie et al. (2002a). 

2.3 Spectroscopy 

We obtained spectra of 2 BD candidates using ISIS on the 
4.2 m William Herschel Telescope (WHT) during 29 Jan 
2001. The seeing was ~1.2 arcseconds and the weather was 
clear. Our instrumental setup comprised the TEK4 CCD 
and the R158R grating on the ISIS red arm, with the mirror 
in place and no order blocking filters, so as to minimize flux 
losses (we expect second order contamination for these red 
objects to be negligible). We used a 1 arcsecond slit, and 
integrated for 1800s. This gave us a wavelength coverage of 
6213-9193A with a dispersion of 2.91 A pixel" 1 . 

We used IRAF to reduce our spectra. The CCD frames 
were bias-subtracted and flat-fielded using standard rou- 
tines, and the spectra extracted and sky-subtracted using 
the apall routine. We calibrated the wavelength scale using 
the CuAr arc lamp exposures, and the flux using observa- 
tions of flux standard stars (available in the IRAF environ- 
ment) on the same night, and with the same instrumental 
setup. 



3 RESULTS & DISCUSSION 
3.1 Pleiades spectra 

Figure 1 shows the spectra of the two Pleiades BD can- 
didates BPL 327 and BPL 45, where they have been nor- 
malised to unity at 7500A. Table 1 shows the spectroscopic 
information extracted from these data. Our adopted spec- 
tral types were derived using the Pseudo-continuum ratio 
PC3 from Martin, Rebolo & Zapatero-Osorio (1996), and 
have uncertainties of ±0.5. For BPL 327 the VO index of 
Kirkpatrick, Henry & Simons (1995) agrees with the PC3 
index. Also, the A-index of Kirkpatrick, Henry & McCarthy 
(1991) implies a slightly earlier spectral type for both sources 
(dM6.5 and dM7.5 for BPL 327 and BPL 45 respectively), 
which is what we expect for young sources that are not fully 
contracted. We also note that the BPL 327 spectra shows 
that the 12 and 13 indices of Martin & Kun (1996) (which 
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Figure 1. The 2 Pleiades BD candidate spectra (normalised to 
unity at 7500A). The noisier spectra of BPL 45 has been binned 
by 2 pixels in this plot, and offset for clarity. 

measure the CaH and the TiO absorption bands respec- 
tively) continue to increase out to ~dM7, before turning over 
and decreasing to lower values by dM8. By comparing the 
I magnitudes and derived spectral types with the spectral 
type- magnitude diagram of Martin et al. (2000), we found 
that both BPL 327 and BPL 45 lie on the Pleiades sequence, 
and appear to be single BDs. 

3.2 The IK CMDs 

Figures 2 and 3 show the K, I-K CMDs (IK CMDs) for the 
Pleiades and Praesepe respectively. Small filled circles are 
higher mass candidate members with 2MASS photometry, 
and proper motion candidates from HSHJ in the Pleiades 
and Praesepe respectively. We have several additional fainter 
Praesepe candidates with 2MASS photometry which are not 
plotted, since their errors are sizable (see Table A. 4). Larger 
filled symbols are lower mass candidates (candidate BDs 
in the Pleiades), where circles have J- and H-band pho- 
tometry (triangles do not). Open circles (small and large) 
look like background sources in these plots. Open diamonds 
and open squares are likely old-disk non-members from the 
K, J-K CMDs (JK CMDs) and the J-H, H-K 2-colour di- 
agrams (JHK diagrams) respectively (see Sections 3.4 and 
3.5). Circled objects look like unresolved binaries in these 
plots (see below). Objects over-plotted with an additional 
square or triangle are likely unresolved binaries from the JK 
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CMDs and JHK diagrams. Typical photometric error bars 
are shown in the bottom left. The location of several magni- 
tude bins are indicated with dotted lines (see Section 3.5). 
Solar metallicity NextGen and Dusty isochrones are over- 
plotted as dashed lines, with mass points shown in units of 
0.001M Q . We also show higher mass (>O.7M0) isochrones 
derived with a larger mixing length parameter (l=1.9H p ), 
since this parameter makes a significant difference at these 
higher masses (Chabrier & Baraffe f 997). In the Pleiades, we 
assume an age of 1 25Myrs, A/=0.06, Ax =0.05, and use the 
re-computed Hipparcos distance modulus of (m-M)o=5.57 
(Makarov 2002), which was derived in a way that reduces 
the propagation of along-scan attitude errors (a source of 
discrepancy with the previous Hipparcos Pleiades distance 
(eg. van Leeuwen 1999). In Praesepe, we show both 0.5 and 
I Gyr isochrones, and assume (m-M)o=6.0 with zero extinc- 
tion. This distance modulus was obtained from a best-fit 
of the NextGen isochrones to the good Praesepe candidates 
with UKIRT photometry from I-K=2.5-3.5. In this colour 
range the 0.5 and lGyr isochrones lie on top of each other. 
This distance agrees well with other values found from main 
sequence fitting (Lynga 1987, Nissen 1988), but is slightly 
lower than the Hipparcos value of 6.24 (Mermilliod et al. 
1997). However, it is possible that the previously mentioned 
attitude errors could affect the Hipparcos distance of Prae- 
sepe. The over-plotted solid lines represent our selected clus- 
ter sequences (constructed both from isochrone data and 
empirical data), with their equal mass binary sequences 
shown as dotted lines (see subsequent discussion). 

We will first discuss the brighter sources with 2MASS 
K-band photometry. For the range K=f0.5-I4 in the 
Pleiades, and K=lf-I4.5 in Praesepe there is a known kink 
in the main sequence, where it becomes steeper. This fea- 
ture is seen over a wide range of wavelengths (Reid & Cruz 
2002), and is thought to result from H2 formation in the 
photosphere, lowering the adiabatic gradient and causing 
convection (Kroupa Tout & Gilmore 1990). The resulting 
lower temperature gradient then leads to increased T e //s, 
with the consequent reduction in the I-K colour, and turn 
over on the CMD. The kink is clearly seen in both cluster 
CMDs. For Praesepe, the model isochrone follows the data 
very well (to within the photometric errors). However, for 
the Pleiades the data appear to be somewhat redder than 
the isochrones from K=f 0.5-13 (Mx~5-7.5). In order to de- 
termine if this effect is an artifact of the different photomet- 
ric filter systems used, we have examined our photometry 
more closely. Photographic magnitudes (used in HHJ) and 
KittPeak magnitudes (used in Pinfield etal 2000) were both 
available for the Pleiades sources in this range, and by trans- 
forming both onto the Cousins system we looked for any 
systematic differences that could account for such an effect. 
However, no such differences were found, and the Pleiades 
sources looked redder than the isochrones whichever I-band 
data was used. The Pleiades and Praesepe sources in this 
range all have 2MASS K-band photometry (converted onto 
the MKO system), and so there should be no systematic dif- 
ference in this band. It thus seems unlikely that this effect 
is due to photometric inconsistencies. 

Assuming this effect is genuine, we should expect to see 
it amongst other similar populations. The open cluster aPer 
has a similar age to the Pleiades (~80Myrs), and Figure 4 
shows an IK CMD of likely aPer members (filled circles) 
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Figure 2. The K, I-K CMD of the Pleiades. Small filled cir- 
cles have 2MASS photometry. Larger filled symbols have UKIRT 
photometry; circles have JHK, triangles have just K. Open circles 
are non members. Circled objects are IK binaries. Objects over- 
plotted with an additional open square or triangle are JK or JHK 
binaries respectively (see text). Dashed lines are the NextGen and 
Dusty models. Solid and dotted lines are the cluster single and 
binary star sequences respectively. 

from Prosser etal (1992) and Barrado y Navascues (2002). 
These have been corrected for distance and reddening (m- 
M=6.23, A/=0.17, A K =0.04; Stauffer etal 1999) for easy 
comparison with the field stars that are also plotted (open 
circles; Reid & Cruz, 2002). NextGen isochrones of 80Myrs 
and I Gyr are indicated by solid lines, and a dotted line in- 
dicates the expected spread due to binarity for oPer. It is 
clear that (despite some scatter) the lGyr model isochrones 
represents the field star data very well. However, the major- 
ity of the cvPer sources from Mk ~5-7 are redder than both 
the single and the binary 80Myr isochrones, and the same 
effect is therefore evident in aPer aswell as the Pleiades. 

So what could be causing this reddening? Variability 
could make a source appear redder since the I- and K- 
band measurements were not made simultaneously. How- 
ever, sources could equally be flaring when the I- or K-band 
measurements were made, and the overall effect should be 
to simply broaden the sequence slightly, not shift it to the 
red. The magnetic activity of such youthful sources may, 
however, provide the explanation for the reddening. Haw- 
ley, Tourtellot & Reid (1999) found that dMe stars in NGC 
2516 (~Pleiades age) lay redward of the fiducual main se- 
quence in an My verses V-I CMD. They also found that 
some molecular bands respond strongly to the presence of 
a chromosphere, and concluded that the V-I colour was af- 
fected by the presence of activity. The I-K colour and the 
V-I colour could be reddened by the same process. It may 
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Figure 3. The K, I-K CMD of Praesepe. Symbols are the same 
as in Figure 2. 
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Figure 4. The K, I-K CMD of Alpha-Per (filled circles) and 
field stars (open circles). Sources arc taken from Prosser (1992), 
Barrado y Navascues (2002), and Reid & Cruz (2002). 



seem surprising that photospheric bands (and the photomet- 
ric colours they influence) should be affected by the presence 
of a chromosphere. However, modeling the correct form of 
the magnetic field that threads the atmospheres of these 
stars may allow us to account for these effects. 

Below the kink the NextGen isochrone fits both cluster 



CMDs well out to I-K~3.5, showing the expected decreasing 
T e //s as degeneracy sets in. For Praesepe, we defined our 
cluster sequence out to I-K~3.5 with the NextGen model 
isochrone, since it fits the data very well in this range. We 
chose to do the same for the Pleiades despite the sources be- 
ing slightly redder than the isochrones for K=10.5-13, since 
we do not expect the reddening to significantly effect the 
mass-magnitude model predictions. 

At redder colours we confine our discussion to the 
Pleiades, since we only have a few good Praesepe candidates 
in this range. It can be seen that there is a fairly obvious gap 
in the Pleiades sequence from K~14. 65-14. 95 (at I-K~3.5) 
where the number of sources is severely reduced. The spec- 
tral type corresponding to the gap is M7 (see Figure 4 of 
Martin et al. 2000), and we refer to this feature as the "M- 
dwarf gap". It is discussed in more detail in Dobbie et al. 
(2002c). However, to summarise, we believe that it results 
from a sharp local drop in the magnitude-mass relation, pos- 
sibly caused by the formation of dust in the atmospheres of 
objects in this T e // regime. We measured the size of the 
gap at 0.3 in K which compares to ~0.5 in I (Dobbie et al. 
2002c). 

Below the Pleiades M-dwarf gap (K~14.95-15.3) the 
single source candidates follow a fairly flat path in the CMD 
out to I-K~4.1. There is then a turn down, where the candi- 
dates drop onto the Dusty isochrone, joining it somewhere 
from I-K~4.1-4.5. We estimated the path followed by the 
flat part of the sequence by averaging the colours and K- 
band magnitudes of two groups of sources; the 3 bluest single 
sources and the 3 reddest single sources from K=14.95-15.3 
(below which the turn down is obvious). And we estimated 
the colour where the sequence joins the Dusty isochrone as 
I-K~4.3. The Dusty isochrone then agrees well with the 5 
faintest Pleiades candidates out to I-K~6, so in this range 
we used the Dusty isochrones for our sequence. 

For Praesepe, we have extended the cluster sequence 
beyond the NextGen isochrones by assuming it has the same 
form as the Pleiades sequence; the same M-dwarf gap size, 
and the same flat portion that turns over and falls onto 
the Dusty isochrone by I-K~4.3. This may not be entirely 
accurate, but should act as a reasonable approximation. 

We have flagged non-members in the Pleiades and Prae- 
sepe using our cluster sequences as a guide, and allowing 
for photometric uncertainties as well as the depth effect 
of the clusters (±0.2 mag for the Pleiades, and ±0.15 mag 
for Praesepe; Pinfield, Jameson & Hodgkin 1998, Holland 
et al. 2000). We have also flagged likely unresolved binary 
members with UKIRT photometry, if they lie significantly 
(>0.25 mag ) above the single star sequence. We duly identi- 
fied 13 Pleiades binaries, and 8 Praesepe binaries in this way. 
We refer to them as IK binaries, and they are overplotted 
with an open circle in the figures. We note that there is a de- 
generacy between mass and binarity where the IK sequence 
is dropping down onto the Dusty isochrone (I-K=3.9-4.5), 
and it is not possible to identify all binaries in this region 
from the IK CMD alone. 

3.3 Model J-K colours 

In the previous section we found that the NextGen 
isochrones fit the IK CMDs very well for I-K~l-3.5, except 
for Pleiades sources from K=10.5-13, which are slightly red- 



© 2003 RAS, MNRAS 000, 1-14 



6 D. J. Pinfield et al. 




Figure 5. The J-K, I-K two-colour diagram of field star sources. 
Filled circles are young disk, and open circles are old disk. Points 
have no age information available. 125Myr and lGyr NextGen 
and Dusty isochrones arc shown as dashed and dotted lines re- 
spectively. 



der than the isochrones. We also found that Dusty isochrones 
fit the Pleiades IK CMD well for I-K=4.3-6. In this section 
we test the J-K model predictions for these ranges, for a 
sample of field stars with kinematics and well measured I-K 
and J-K colours. 

Figure 5 shows J-K plotted against I-K for sources taken 
from Leggett (1992), Leggett, Allard & Hauschildt (1998) 
and Dahn et al. (2002). All sources classified as young or old 
disk (from their kinematics) in the Leggett papers are shown 
as filled or open circles respectively. In addition, we have 
assigned 4 L dwarfs as young (<0.5Gyrs), based on lithium 
detections (Martin et al. 1999). These are also shown as 
filled circles. The uncertainties are generally ±0.05 in colour. 
We have also overplotted the NextGen and Dusty model 
isochrones for 120Myrs (dashed lines) and lGyrs (dotted 
lines) . 

Comparison of the two populations shows that for bluer 
sources (I-K<3.5, J-K<1), the old disk stars have a bigger 
spread in J-K and on average are slightly bluer in J-K than 
their young disk counterparts. The two NextGen isochrones 
essentially lie on top of each other in this range, and fit the 
data very well. For redder sources (I-K=4.3-6, J-K=l.l-2), 
the Dusty isochrones lie on top of each other, and are in 
excellent agreement with the young disk data. In between 
these two ranges the isochrones are all too blue in J-K. The 
NextGen models are too blue because they do not account 
for dust, and the Dusty models are too blue because they 
use a less sophisticated treatment of H2O and TiO molecular 
opacities. This is all consistent with the gradual change- 
over from NextGen to Dusty isochrones that was seen in the 
IK CMDs, and suggests that these models should provide 
accurate NIR isochrones in the blue and red photometric 
ranges. 



Figure 6. The K, J-K CMD of the Pleiades. Symbols are the 
same as in Figure 2. 
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Figure 7. The K, J-K CMD of Praesepe. Symbols are the same 
as in Figure 2. 



3.4 The JK CMDs 

Figures 6 and 7 show the JK CMDs for the Pleiades and 
Praesepe respectively. Only good candidates from the IK 
CMDs are shown. The symbols mean the same as in Figure 
2. The location of several magnitude bins are indicated with 
dotted lines (see Section 3.5). Isochrones are overplotted as 
dashed lines, with mass points shown in units of O.OOIMq. 
Typical photometric errors for UKIRT photometry are in- 
dicated. 

We constructed cluster sequences for the JK CMDs 
based on our analysis in Sections 3.2 and 3.3. We used 
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NextGen isochrones for I-K<3.5 (mass<0.07 and O.O9M 
for the Pleiades and Praesepe respectively), and Dusty 
isochrones for J-K=l.l-2 (I-K=4.3-6). In between the 
NextGen and Dusty isochrones we added two additional se- 
quence points, based on our Pleiades sources; one to bridge 
the M dwarf gap, and the other to represent the turn down 
point, where the sequences begins to fall onto the Dusty 
isochrone (see Section 3.2). An average of the J-K colours 
of the Pleiades sources above and below the gap shows that 
there is essentially no significant J-K colour change across 
the gap (AJ-K=0.01). Averaging of the 3 single Pleiades 
sources in the turn down region gave a turn down colour of 
J-K=1.03. The flat region seen in the IK CMDs is thus not 
very apparent in the JK CMDs. For Praesepe we assumed 
that the JK colour and magnitude changes are the same as 
in the Pleiades (as we did previously for the IK CMDs). The 
cluster sequences we constructed are shown as solid lines in 
the JK CMDs, and a dotted line indicates the spread we 
expect in the sequence from unresolved binarity. 

The Net Gen part of the sequence fits the observations 
very well in both clusters, and there is certainly no evidence 
for the Pleiades sources being redder than the isochrones 
from K=10.5-13, as was seen in the IK CMDs. The J-K 
colour is clearly not affected in the way that the I-K colour 
is. Below the M-dwarf gap, it can be seen that there is 
a spread in the Pleiades sources, where several are signif- 
icantly redder than expected for single objects. A spread 
is also seen for a small number of Praesepe sources. Such 
a spread could result from several possibilities. Transient 
surface features such as spots could be responsible. The 2 
most obvious red sources (Teide 1 and Calar 3; BPL 137 and 
235) have been tested for variability. Bailer- Jones & Mundt 
(2001) found no evidence of I-band variability in Teide 1 
over a 100 hour time scale, and only marginal evidence (just 
above the noise) for sporadic variability (AI<0.03) in Calar 

3 over a 29 hour time scale. Our photometry of Teide 1 
comprises 3 well separated (>a year) J- and H-band mea- 
surements which are consistent to within their errors. For 
Calar 3 we have 1 K-band measurement which is consistent 
with the measurement presented by Martin et al. (2000), 
and 2 H-band measurements which agree to within their 
errors. We also have 3 J-band measurements which show 
marginal evidence (at the 2a level) of variability. The level 
of this variation is not in itself sufficient to account for the 
redness of Calar 3. Photometric errors will also cause some 
spread in the sources. However, the general consistency (at 
the 1-a level) between our multiple measurements makes it 
unlikely that these alone are responsible. The third option, 
and the one we chose to consider in more detail, is that addi- 
tional unresolved binaries (not identified in the IK CMD) are 
present. A substantially redder unresolved companion could 
redden a source sufficiently to account for the spread. These 
binaries would either be background sources, or cluster bi- 
naries that come from the turn down region of the IK CMD, 
where binarity is not always apparent. We have thus flagged 

4 additional candidate unresolved binaries in the Pleiades, 
and 2 in Praesepe. We refer to these as JK binaries, and 
they are overplotted with an additional open square in the 
figures. We note that the bluest of the Pleiades JK binaries 
is a known background binary (with out lithium; Martin et 
al. 2000), resolved by the HST (Martin et al. 1998). The 
Pleiades binaries will be discussed in more detail in Section 



3.7. Of the 2 Praesepe JK binaries, the bluer one (RIZ Pr 23) 
is located just above where we expect the missing M dwarf 
gap to be, and could be a binary from below the gap. The 
redder one (RIZ Pr 23) comes from the turn down region of 
the IK CMD. 

We have flagged non-members in the Pleiades and Prae- 
sepe using our cluster sequences as a guide. Objects that 
look too blue are presumably background stars from an old- 
disk population (see Section 3.3), and we flagged them as 
such if their J-K colours were at least 2-a blueward of the 
cluster sequence. 

3.5 The J-H, H-K 2-colour diagrams 

Figures 8 and 9 show the J-H, H-K 2-colour diagrams (JHK 
diagrams) for the Pleiades and Praesepe respectively. Only 
good candidates from the IK CMDs are shown. The source 
symbols mean the same as in Figure 2. The positions of 4 
field L-dwarfs (from Dahn et al. 2002) are also shown as as- 
terisks. The photometric error bar shown in the lower right 
of the plots are typical of the photometry of the BPL sources. 
The fainter sources have slightly larger errors. Solar metal- 
licity model atmosphere predictions (used for the NextGen 
and Dusty models) are overplotted as dashed and solid 
lines respectively, with labels indicating surface gravities. 
For the Pleiades we plot NextGen models for T e // >2000 
and 3000K, and Dusty models for T e// <2000K. For Prae- 
sepe, we plot NextGen models for *T e ff >2500 and 3000K, 
and Dusty models for T eff <2500K. 

For the highest mass stars, it can be seen that the 
NextGen models are ~0.1 too blue in H-K up to the 2-colour 
turn-over, where convection begins (~4000K). The data 
turns over at H-K=0.27, J-H=0.57, and log g for such objects 
should be ~4.5. The theoretical J-H colour of the turn over 
is therefore ~0.1 too red. These discrepancies are thought to 
be caused by inaccuracies in the modeling of water vapour 
opacities in the H-band (with its overall low opacity) as well 
as an underestimation of the mixing length parameter. Af- 
ter the turn-over, the 2MASS photometry follows a short 
path in the JHK diagram, reaching J-H=0.52, H-K=0.38 by 
~3200K, in reasonable agreement with the log g—5 NextGen 
model. In order to look for any non-evolutionary differ- 
ences between the colours of Pleiades and Praesepe stellar 
members, we selected Pleiades members from K=11.5-12, 
and Praesepe members from K=12-12.5. These two sets of 
sources should be fully contracted (according to the models), 
and should represent the T e // ~3500-3600K range (since 
we have corrected for distance). We see no significant differ- 
ence in the colours, with J-H=0.53±0.03 and 0.54±0.03, and 
H-K=0.35±0.03 and 0.33±0.04 for the Pleiades and Prae- 
sepe respectively. Praesepe is known to be slightly metal 
rich compared to the Pleiades, with [Fe/H]=0.04±0.04-0.14 
(Friel & Boesgaard 1992; Reglero & Fabregat 1991). How- 
ever, this does not appear to significantly affect the stellar 
NIR colours, consistent with the model predictions in this 
T eff range. 

In order to show what is happening at cooler T e //s, we 
grouped the sources into a number of K magnitude bins. 
Above the gap we selected Pleiades bins of K=13. 6-14.1 
and K=14.1-14.65, and Praesepe bins of K=14.5-15.3 and 
K=15.3-15.85. We labeled these bins A and B respectively, 
and they should correspond to approximately the same T e / / 
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Figure 8. The J-H, H-K 2-colour diagram of the Pleiades. Sym- 
bols are the same as in Figure 2. Individual sources are labeled 
in italic (i?=Roque). 
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Figure 9. The J-H, H-K 2-colour diagram of Praesepe. Symbols 
are the same as in Figure 2. 

ranges in each cluster (-3150-3000K for A, and 3000-2800K 
for B). Below the gap, we considered 3 additional Pleiades 
bins. The K=14.9-15.13 range (bin C) contains objects 
on the flat section of the IK CMD (see Section 3.2). The 
K=15.13-15.5 range (bin D) contains objects at the top of 
the turn down, where the sequence begins to drop towards 
the Dusty isochrone in the IK CMD. Finally, the K=15.5- 
16.1 range (bin E) contains objects from near the base of 
the turn down region, where the IK sequence is joining onto 
the Dusty isochrone. These bins are indicated on the IK and 
JK CMDs with dotted lines, and in the JHK diagrams by 



boxes drawn with dotted lines. In the Pleiades, bin A con- 
tains 2MASS sources, and we defined the box dimensions 
using the standard deviation of the average bin colours. For 
the remaining bins we defined the box dimensions by at- 
tempting to include all non-binary member sources within 
the box dimensions. 

When we compare the colour-colour location of bins A 
and B in each of the clusters, it can be seen that they differ. 
In the Pleiades, bins A and B have similar colours, and lie 
very close to the NextGen model predictions for their T e ff 
range. However, in Praesepe these bins have J-H colours 
that are ~0.1 bluer than this, and bin B has an H-K colour 
~0.1 redder. The surface gravities expected for these bins 
are logg=5 for the Pleiades, and log g=5. 15-5.25 for Prae- 
sepe, and these differing surface gravities may provide the 
explanation for the differing colours. We do not expect any 
atmospheric dust formation for objects in these bins since 
their T e //s are too high (see Tsuji 2002 figure 7, and Jones 
& Tsuji 1997), so gravity sensitive dust formation cannot 
provide the explanation. However, molecular opacities (such 
as H2O) could be more sensitive to gravity (or pressure) 
than the NextGen models predict and future generations of 
models with updated EOS and opacity data may account 
for these colour changes. 

Moving to fainter cooler objects in Praesepe, it can be 
seen that the faint binary candidate below bin B lies on 
the dusty isochrone, consistent with dust formation begin- 
ning in Praesepe at around T e ff ~2500K. However, in the 
Pleiades we are faced with a rather different story. Mov- 
ing across the gap from bin B to C produces no significant 
change in colour, and between bins C and D there is then 
an increase of ~0.1 in J-H. These changes are all consis- 
tent with the NextGen isochrone. Pleiades sources eventu- 
ally move onto the dusty isochrone between bins D and E, 
and the 2 faintest Pleiades sources (Roque 25 and 30) both 
lie on the Dusty isochrones. A possible explanation for these 
near infrared colours is that dust does not begin to form in 
Pleiades sources until T e / / ~2000K (between bins D and E). 
This might occur because in lower gravity atmospheres the 
lower pressures could suppress dust formation until lower 
T e //s. However, one would also expect that lower gravity 
atmospheres would be more extended, and have lower outer 
temperatures, which would have the opposite effect. 

Assuming that the differences in colour between 
Pleiades and Praesepe sources result from differing surface 
gravities, then one might expect even greater differences for 
sources found in star forming regions. This would have im- 
plications for the derreddening process that is commonly 
used for such sources (eg. Lucas & Roche 2000). 

We have flagged likely non-members using the location 
of bins A to E and the Dusty model atmosphere predic- 
tions (for the faintest sources) as a guide. Old disk sources 
appeared too blue in J-H, and several non-stellar sources 
were also identified amongst our data by their extremely red 
colours. One object in Figure 7 (Roque 33) appears to be too 
blue in J-H. This object is a proper motion member how- 
ever (see Table A.l), and we note that its photometric errors 
are ~0.2 in colour making them not inconsistent with mem- 
bership. We also used the Pleiades JHK diagram to identify 
some additional unresolved binaries. In the turn down region 
of both the IK and the JK CMDs (where the sequence is 
dropping steeply down towards the Dusty isochrones) , there 
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remains some degeneracy between mass and binarity. How- 
ever, we were able to break this degeneracy using our JHK 
diagram. A single source with K=15. 13-15.5 should lie in 
bin D. However, 3 sources in this magnitude range lie in (or 
very near to) bin E. These are presumably near equal mass 
binaries, where both components of each binary come from 
bin E (we refer to them as JHK binaries). These binaries are 
overplotted with an additional open triangle in the figures. 

3.6 Summary of membership criteria 

Using the IK and JK CMDs, as well as the JHK diagrams, 
we have assigned photometric membership criterea to our 
Pleiades and Praesepe samples (see Appendices). In the 
Pleiades, we flagged 5 of the 30 BPL BD candidates as non- 
members. The remaining 25 BPL candidates, and the 9 addi- 
tional Pleiades candidates from the literature look like clus- 
ter BDs. Examination of Table Al shows that our criterea 
for these sources are consistent with existing spectroscopic 
and astrometric measurements, where they exist. It is par- 
ticularly encouraging that we were able to flag CFHT-PL- 
18 (BPL 283) as a background binary from its photometry 
alone. This source was a suspected BD binary, and was only 
confirmed as a background source when Martin et al. (2000) 
failed to detect lithium in its spectrum. BPL 45 has been 
flagged as the same type of object based on its measured 
J- and K-band photometry alone. Of the Pleiades low-mass 
star candidates, we flag 14 non-members out of the 194 can- 
didates. This represents only 7% spurious sources in the 
sample. In Praesepe, we identify higher levels of contami- 
nation for our faintest sources. We flag 20 non-members out 
of the 26 RlZ-Pr candidates, and 37 non-members out of 
the 89 IZ-Pr candidates. However, for the brighter HSHJ 
sources we only flag 21 non-members out of 458. The posi- 
tion of Roque-Pr 1 in the IK and JK CMDs is consistent 
with membership. However, if it is a member, it must be an 
unresolved BD binary. 



3.7 BD binary mass-ratios 

Figure 10 shows two regions of the Pleiades IK and JK 
CMDs, focusing on the candidate BDs. The symbols mean 
the same as in Figure 2. NextGen and Dusty isochrones are 
shown with dashed lines, and our cluster sequences (see Sec- 
tions 3.2 and 3.4) are shown as solid lines. We also show the 
effects of unresolved binarity, with the use of binary tracks 
(shown as dotted lines). To construct these binary tracks 
we started with the single star point which is to form the 
primary. Then, by combining the magnitude of this primary 
with that of the faintest point on the sequence (lowest mass) , 
we obtained a new binary point which we joined to our pri- 
mary point by a dotted line. We then repeated this process 
for higher mass secondaries, joining each new point to its 
predecessor until we reach an equal mass binary, 0.75 mag- 
nitudes above the original primary. We also show 5 mass 
points on the CMDs (in units of 0.001M©). The O.O7M 
point just above the M dwarf gap is the NextGen model 
prediction for a 125Myr isochrone. Since we believe that the 
M dwarf gap results from a sharp drop in the luminosity- 
mass relation, we assumed that the mass just below the M 
dwarf gap is very similar to the mass just above it. Hence, we 



have a O.O7M0 point just below the M dwarf gap also. The 
0.06M© mass point was estimated from the Pleiades mass-I 
relation of Dobbie et al. (2002c), which accounts for the M 
dwarf gap, and was derived using the Pleiades I-band BD 
LF, and an assumed power-law MF. We note that this agrees 
with the NextGen and Dusty model J-band predictions, if 
they are offset by the J magnitude width of the M dwarf gap. 
The final two mass points are Dusty model predictions for 
a 125Myr isochrone. The O.O4M0 mass point should be rea- 
sonably representative of the point where the sequence joins 
the Dusty isochrone. We also note that Dobbie et al. (2002c) 
suggest that model predictions for the masses of Pleiades 
BDs from 1=18.5-19.5 are too low (by up to 0.01M©), which 
may have some bearing on the Dusty mass estimates. 

When analysing these CMDs, we must consider the 
depth effect of the cluster. Sources near the front and back 
of the cluster will appear slightly brighter and fainter respec- 
tively, than they would in the cluster centre. The tidal ra- 
dius of the Pleiades is 13.1pc (Pinfield, Jameson & Hodgkin 
1998) giving a full depth of ±0.2 ma9 . However, the majority 
of sources will be found closer to the cluster centre. Jame- 
son et al. (2002) fit a King profile (King 1962) to the spatial 
distribution of Pleiades BDs, and derived a core radius of 
5pc. We would expect ~80% of members to lie within this 
core radius, giving a more likely depth of ±0.1 ma9 . In order 
to estimate binary mass-ratios (q) for our unresolved binary 
candidates, we compared their positions to the binary tracks 
in both CMDs, and estimated q ranges consistent with the 
photometric errors and the depth effect of the Pleiades. We 
did this for a depth effect of ±0.1 and ±0.2. We also esti- 
mated q ranges on the basis that the Dusty isochrone masses 
may be higher by 0.01 M©. And finally, we made a small al- 
lowance for the known (or suspected) variables (PPL1 and 
Calar 3). Our results are summarised in Table 2. Column 
3 indicates if the source is an IK, a JK or a JHK binary 
(identified from the IK CMD, the JK CMD of the JHK dia- 
gram) . Column 4 gives the estimated q range from the mass 
points in Figure 9 and a ±0.1 depth effect. Column 5 gives 
a larger q range, estimated with the full ±0.2 depth effect, 
and accounting for the uncertainty in the Dusty isochrone 
mass points (assuming they may be 0.01M© too low). 

Where a source appeared in both CMDs, we found q 
ranges that were consistent (to within photometric errors) 
for all but 2 candidates; CFHT-PL-18 (BPL 283) and BPL 
45 appear red in the JK CMD, but do not appear sufficiently 
red in the IK CMD. The component sources of these binaries 
would therefore be too blue in I-K to be cluster members. 
This is confirmed for CFHT-PL-18, which is a known back- 
ground binary (see Section 3.4). We have also flagged BPL 
45 as a probable non-member, based on its measured pho- 
tometry alone. CFHT-PL-12 (BPL 294) and INT-PL-IZ-42 
could have primaries that are either just above or just below 
the M dwarf gap. However, CFHT-PL-12 would need to be 
near the front of the cluster if its primary was from below 
the gap. Roque 4 lies quite high in the JK CMD suggesting 
that it may be a foreground object. However, it could still 
be a binary BD if it is at the front of the cluster, and may 
alternatively be a BD triple system. We were unable to es- 
timate good q ranges for the 2 faintest binaries, since we do 
not know the form of the cluster sequence for masses below 
0.03M©. 

It is clear that the majority of the unresolved binary 
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Table 2. Pleiades BD binary mass-ratios (q). The qi ranges were estimated assuming a depth effect of ±0.1 ma9 , and masses from Figure 
9. q2 ranges were estimated assuming a depth effect of ±0.2 ma9 , and accounting for the fact that the Dusty isochrone masses could be 
underestimated by -~0.01Mq. 



Binaries Other names Binary type qi q2 Notes 



BPL 45 




JK 






Probable BG source 


BPL 66 


Roquc 4 


IK 




~1 


Near front of clustcr?/Triple BD? 


BPL 79 


Roquc 13 


IK 


0.7-1 


0.65-1 




BPL 108 


Roque 14 


IK 


0.8-1 


0.75-1 




BPL 137 


Teide 1 


JK 


0.55-0.65 


0.5-0.9 




BPL 142 


Roquc 17 


IK 


0.7-1 


0.65-1 




BPL 235 


Calar 3 


JK 


~0.5 


0.5-0.7 




BPL 249 




JHK 


~1 






BPL 283 


CFHT-PL-18 


JK 






Known BG source 


BPL 294 


CFHT-PL-12 


IK 


0.6-0.9 


0.4-1 


Primary probably above gap 


BPL 303 


CFHT-PL-25 


JHK 


<0.75 


<1 




PPL 1 


Roquc 15 


IK 


0.8-1 


0.75-1 




PPL 15 


IPMBD 23, NPL35 


IK 


0.7-1 


0.65-1 


Known binary with q=0.85 


CFHT-PL-16 




IK 


0.8-1 


0.75-1 




CFHT-PL-23 




JHK 


~1 






INT-PL-IZ-20 




IK 




<0.8 




INT-PL-IZ-25 




IK 


~0.8 






INT-PL-IZ-42 




IK 


0.4-0.6 or ~1 


0-0.7 or ~1 


Primary above or below gap? 


INT-PL-IZ-43 




IK 


0.4-0.65 


0-0.8 




INT-PL-IZ-76 




IK 


0.9-1 


0.8-1 





candidates have fairly high q values. Indeed, 15 of the 18 
probable binary BDs appear to have q>0.7. However, the 
photometry of 3 sources is inconsistent with them being ei- 
ther single or q>0.7 binaries. One of these is without J- or 
H-band measurements, and has not been confirmed as a clus- 
ter member (by either a proper motion measurement or the 
lithium test). The other 2 are Teide 1 and Calar 3, which are 
bona fide cluster BDs. We have previously discussed (Sec- 
tion 3.4) other possible causes for the red J-K colours seen 
in these sources, and concluded that binarity was the most 
likely. We note in addition that the position of these sources 
in the JHK diagram is also consistent with binarity, since 



they appear partway between bins D/E and the location of 
Roque 25. These few objects therefore suggest that BD bi- 
naries are not all high q systems. They could consist of two 
populations; the majority in high q systems, and a smaller 
fraction in low q systems. Theoretical simulations are not 
yet at a level to make such predictions. However, the hy- 
drodynamical simulation of Bate, Bonnell & Bromm (2002) 
shows BDs forming from the collapse and fragmentation of 
a turbulent molecular cloud. They determined that ^75% 
of BDs form by the fragmentation of gravitationally unsta- 
ble circum-stellar disks, with subsequent ejection removing 
the BD from the region of dense gas. The remaining 25% 
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Table 3. Binary fractions for the Praesepe and Pleiades samples. 
Spectral types were estimated using Leggett (1992). 



Praesepe 



I-K Mass/SpT q-sensitivity Binary fraction 



1-1.9 


1.0-0.6M Q 
G/K stars 


0.5-1 


17jl|% 


1.9-2.5 


0.6 — O.2M 
M0-M5 


0.35-1 


41±5% 


1.9-2.2 


0.6 — O.35M 
M0-M3 


0.3-1 


31± 7 6 % 


2.2-2.5 


~0.35 — O.2M 
M3-M5 


0.4-1 


44±6% 


3-3.6 


~0.11 — O.O9M 
M6-M7 


0.65-1 


47±\l% 


Pleiades 


1-1.9 


1.0-0.6M Q 
G/K stars 


0.5-1 




2.6-3 


O.15-O.O9M 
M5-M6 


~0.35-l 


36±5% 


2.6-2.8 


0.15-O.HMq 


0.35-1 


27±*% 




M5 




2.8-3 


O.ll-O.O9M 
M6 


0.4-1 


51 ±8% 


3.3-4.3 


O.O7-O.O6M 
M7-M8 


0.5-1 


50±\l% 



were found to form from the fragmentation of dynamically 
collapsing filamentary gas, where the forming BDs fell into, 
and were rapidly ejected from unstable multiple systems (see 
also Reipurth & Clarke 2001). However, this simulation only 
produced 1 BD binary, but lacked the spatial resolution to 
predict very close systems. Future simulations with higher 
spatial resolution may reveal a link between the two differ- 
ent formation mechanisms and the resulting q distribution 
of BD binaries. 

3.8 Binary fractions 

We determined binary fractions (BF; defined as the number 
of binary systems in some photometric range divided by the 
total number of systems in that range) using the IK CMDs 
in general, for a number of photometric ranges. The results 
are given in Table 3, where the uncertainties are determined 
assuming binomial statistics (see Burgasser et al. 2002). The 
mass ranges were estimated using our IK cluster sequences 
(see Section 3.2). The stellar q ranges were estimated by sim- 
ulating the effect of unresolved binarity in the IK CMDs (as 
was done in Section 3.7 for the Pleiades BDs), and finding 
the q value when the binary track crosses midway between 
the single star sequence and the binary sequence 0.75 mag- 
nitudes above it. 

Starting in Praesepe, the first range we considered was 
from I-K= 1—1.9 (G and K stars). The single stars all bunch 
around the predicted sequence, and we flagged binaries by 
their proximity to the binary star sequence. The next range 
we chose was from I-K=l. 9-2.5 (where the upper limit rep- 
resents the completeness of 2MASS and the HSHJ survey). 
The sources lie in a fairly straight path down the CMD in 
this range, and we therefore chose to separate binaries and 



single stars using a straight line. We defined this separation 
line such that the average spread of members above and 
below it was the same (as one would expect for a binary 
induced spread). And by doing this for two sub-ranges (I- 
K=1.9— 2.2 and 2.2-2.5) we ensured that the gradient of the 
line was correct. We would expect a spread of ±0.375 to rep- 
resent the full spread induced from binarity, and we obtained 
spreads of ±0.36 and ±0.43 for the two ranges respectively. 
The slightly larger value for the redder range is accounted for 
by larger photometric errors. This approach has the added 
benefit that it is independent of any assumed distance mod- 
ulus. Table 3 gives the binary fraction for both sub-ranges, 
as well as for the full range. For the lowest-mass Praesepe 
stars we chose to use I-K=3-3.6 (where we just have UKIRT 
photometry). This avoids the CMD overlap region between 
the photographic and CCD surveys, where different survey 
areas and photometric errors complicate matters. Binaries 
in this range were flagged in Section 3.2. Unlike the higher 
mass ranges, these candidates do not have proper motions. 
However, using the disk luminosity function of Kirkpatrick 
et al. (1994), we expect no more than ~1 and ~0.5 field 
star to contaminate the single and binary star sequences 
respectively. Such numbers will not significantly affect our 
statistics. 

For the Pleiades, the first range we considered was from 
I-K=2.6-3. We took the same approach as for Praesepe, and 
divided this range in two (I-K=2.6-2.8 and 2.8-3), defining 
a selection line such that the scatter above and below (in 
each region) was identical. The spread about the selection 
line was ±0.57 for each of the sub-regions. This is somewhat 
more than the expected ±0.375, but the additional spread 
will be predominantly due to photometric errors and clus- 
ter depth. Finally, we considered a BD range. We counted 
single BDs below the M dwarf gap that have K<15.5, since 
our sample is spatially complete in this range. The mass 
range for the single stars is ~0. 07-0. 05 M , where we have 
interpolated between the 0.06 and the O.O4M points in Fig- 
ure 9. We then counted binary BDs that could have pri- 
maries amongst these single sources, including objects from 
K=14.2-15.2. We note that we ignored CFHT-PL-12 (BPL 
294) since it is likely that its primary comes from above the 
gap. In this way, we selected 10 single BDs, and 10 binary 
BDs. We estimated the q sensitivity range for our BD bina- 
ries using Table 2. 

We have already discussed how the cluster depth af- 
fects BD binaries. However, it will also be a source of error 
for our stellar BFs. Objects whose position on the CMD is 
within the cluster depth (eg. 0.2 and 0.15 magnitudes for 
the Pleiades and Praesepe respectively) of the binary selec- 
tion line could potentially spill over this line if they were 
near the front or the back of the cluster. This represents 
53% of Pleiades sources, and 40% of Praesepe sources. How- 
ever, the cluster spatial distributions (Pinfield, Jameson & 
Hodgkin 1998; Holland et al. 2000) indicate that of these 
stars, we would only expect ~5-15% and ~5-10% of Pleiades 
and Praesepe stars respectively to actually spill across the 
line. Further more, an error will only occur if the BF itself 
is different from 50% (for BF=50% misidentified binaries 
and single stars will cancel out). Combining these together, 
we would expect no more than ~1% and 2% BF errors in 
Praesepe and the Pleiades respectively, due to depth effects. 

Another source of error comes from variability, resulting 
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Figure 11. The binary fractions of the Pleiades (triangles) and 
Praesepe (Circles). Three points from the literature arc shown 
as open diamonds. At O.IMq the BF error bars of two points 
overlap, so we plot the Praesepe error bars as a dotted line for 
clarity, q ranges are indicated next to each point, and an A range 
for the Burgasser and Close points. 



from chromospheric activity such as spots or flares pushing 
some single stars intermittently above the binary selection 
line. Lockwood et al. (1984) found that ~40% of Hyades 
F G and K dwarfs varied at a level of ~l-2%. For later 
types, Stauffer et al. (1991) reported H a emission in ~60% 
of late K-mid M Hyades members, aswell as in the vast ma- 
jority of G-mid M Pleiades members, and Scholz, Eisloffel 
& Mundt (2000) found photometric variations among young 
M-dwarfs similar to those of solar type stars. Therefore, al- 
though many of our members may be variable, we expect 
this variability to be generally low (a few %). More violent 
flares are occasionally seen in late M dwarfs (Leibert et al. 
1999), but are rare. We estimate that our stellar BF values 
should not be increased by more than ~5% as a result of 
variability. 

The BF values are plotted in Figure 11, as a func- 
tion of mass. Filled triangles and circles represent Pleiades 
and Praesepe data respectively. Additional points from the 
literature are shown as open diamonds. The highest mass 
point agrees well with the field star results of Duquennoy 
& Mayor (1991), which increases to 70% for q=0.1-l. From 
0.6-0. 35M we find a BF of ~30%. The BF of NGC 2516 
agrees with this value for ~0.6M Q , but drops to ~15% for 
O.5-O.35M (Jeffries, Thurston & Hambly 2001). This dif- 
ference could be explained by the fact that the NGC 2516 
CMD was only sensitive to q=0.6-l binaries, whereas our 
CMD is sensitive to q=0.3-l. If the clusters have intrinsi- 
cally similar binary populations, then this implies that there 
is a BF of 15% from q=0.3-0.6 as well as from q=0.6-l, 
suggesting a flat BF-q distribution in this mass range. This 
represents a departure from the BF-q distribution of solar 
type field stars, which rises towards lower q. If the BF-q 
distribution were like that of solar type stars we would ex- 



pect a BF of 30% for q=0.3-0.6, and thus 45% for q=0.3-l. 
Our result is 1.5a lower than this value. A flat binary-q 
distribution would give a total BF of ^45% in this mass 
range. Next, with the exception of the Pleiades point at 
0.13Mq, there is a general trend of increasing BF towards 
lower mass. The q ranges of these points are fairly similar, 
and so should not significantly bias this trend. Moving into 
the sub-stellar regime, it can be seen that the BD binary 
fraction is around 50%. Even if we are extremely conserva- 
tive, and only select the most obvious BD binaries from the 
IK CMD (assuming that JK and JHK binaries are in fact 
single), we find a BF of 321"%. Figure 11 also shows the 
T dwarf BF resolved by the HST (Burgasser et al. 2003), 
and the late M/early L dwarf BF resolved by Hokupa'a on 
Gemini (Close et al. 2003). Mass range estimates were taken 
from Close et al. (2003) , and estimated for the T dwarfs as- 
suming an age range of l-5Gyrs. Both these BFs (9l4 5 % 
and 15±7% respectively) are significantly lower than our 
BF for Pleiades BDs. However, these observations were only 
sensitive to binaries with separations >1AU. It is possible 
that the total BF of these populations is identical to that 
of Pleiades BDs if 70-80% of BD binaries have separations 
<1AU. It may be that many BD binaries are like PPL15 
(Basri & Martin, 1999) which has a seperation of 0.03AU. 
We also note that the N-body simulations of Kroupa (2001), 
and Kroupa, Aarseth & Hurley (2001) suggest that ^50% 
of cluster BD binaries with log P /days— 1-6 (seperations of 
~0.04-90AU) would be disrupted during the first few Myrs 
of proto-cluster dynamical evolution. This would require an 
even higher initial BD BF unless most BD binaries are ex- 
tremely tight. However, the hydrodynamical simulation of 
Bate, Bonnell & Bromm (2002) predicts an initial BF of 
~5%, since only 1 BD binary was produced in their simu- 
lation among 20 single BDs. But the spatial resolution of 
this simulation was not high enough to model the forma- 
tion of extremely close BD binaries. Such binaries may form 
during a secondary collapse phase of individual pressure sup- 
ported fragments (Bonnell & Bate 1994). If this is not the 
case (Bate 1998), then post fragmentation processes (such 
as dynamical interactions) might be responsible for creating 
extremely tight BD binaries (by hardening) . 



4 FUTURE WORK 

Proper motion measurements have already been made for 
15 of the BD candidates discussed here. Work is currently 
underway by some of us to measure proper motions of the 
remaining Pleiades BDs. It will be particularly useful to 
confirm all possible low-q BD binaries as Pleiades mem- 
bers. Measuring proper motions of the Praesepe candidates 
is also important, and is in progress (Chappelle, Pinfield & 
Steele 2002). Such measurements will allow us to define the 
Praesepe cluster sequence with more confidence. Also, the 
faintest Praesepe candidates could be BDs, and would be 
the first identified in an intermediate age cluster. 

Measuring NIR spectra is a natural next step if we are to 
understand the nature of the cluster sequences. These spec- 
tra could be compared to the predictions of both NextGen 
and Dusty model atmospheres, to determine when dust is 
beginning to form in the atmospheres of Pleiades and Prae- 
sepe sources, and what effects it has in the different surface 
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gravities. Furthermore, the best way to improve on our anal- 
ysis of binarity is from NIR spectra. We could combine the 
spectra of likely single sources, and compare the results with 
spectra of suspected unresolved binaries, to confirm or not, 
if these sources really are binary in nature. 

The nature of the M-dwarf gap could be further in- 
vestigated by measuring longer wavelength photometry (eg. 
L-band) of sources on either side of the M dwarf gap, to look 
for evidence (or lack of it) of dust formation. 

Deeper and larger scale surveys of the clusters will pro- 
vide much larger samples of members to work with. The 
16 square degree survey recently obtained as part of the 
WFS with the INT and WFC should identify several hun- 
dred Pleiades BDs. Also, the WFC survey of Praesepe made 
by Chappelle, Pinfield & Steele (2002) should identify large 
numbers of very low-mass stars and possibly BDs. With 
these larger samples we will be able to more accurately mea- 
sure the M dwarf gap in both clusters, as well as obtain 
a better statistical measure of the BD binary fraction and 
binary-q distribution. 
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Table Al. Photometry (c3-6), membership criterea (c7-9,cll-13) and binarity (clO) of Pleiades BD candidates. Membership is from the 
IK CMD, the JK CMD and the JHK diagram. Refs; (a) Martin, Rebolo & Zapatero-Osorio (1996), (b) Rebolo et al. (1996), (c) Zapatero 
Osorio et al. (1997), (d) Martin ct al. (1998) (c) Stauffcr ct al. (1998), (f) Martin et al. (2000), (g) Pinficld ct al. (2000), (h) Moraux, 
Bouvicr & Stauffcr (2001), (i) Dobbie et al. (2002), (j) Dobbie (2002) private comm. A * indicates a non-member by virtue of binarity. 
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N e 








15.40±0.04 


14.84±0.02 


















BPL 235 


Calar 3 


18.91 9 


16.24±0.03 9 


15.53±0.03 9 


14.91±0.03 9 


y 


y 


y 


b 




yh 


y" 




CFHT-P1-21 




16.09±0.04 


15.49±0.05 


















BPL 240 




18. 68 9 


15.99±0.03 


15.44±0.05 


15.11±0.03 9 


y 


y 


y 










BPL 249 




20. 02 9 


16.66±0.03 9 


16.13±0.03 9 


15.46±0.03 9 


y 


y 


y 


b 








BPL 254 


Teide 2 


17.81 9 


15.53±0.04 


14.94±0.03 


14.57±0.06 9 


y 


y 


y 




fl 
yd 


yh 


y d,e 




CFHT-P1-13 




15.53i0.08 1 


14.97±0.07* 


14.54±0.06 l 
















BPL 283 


CFHT-P1-18 


18. 58 9 


16.04±0.03 9 


15.43±0.02 9 


14.94±0.06 9 


y 


N* 


y 


b 


f 

y' 


N' 1 


AT f 


BPL 294 


f 1 1 — > x x r~r\ x~*1 -1 c\ 

CFHT-P1-12 


17.85 9 


15.16±0.03 


14.55±0.02 


14.20±0.02 9 


y 


y 


y 


b 


y e 


h 

y 


y e 


BPL 303 


CI HT-P1-25 


19. 69 9 


16.63±0.04 


16.01±0.02 


15.47±0.06 9 


y 


y 


y 


b 




yh 










16.64±0.15 8 


16.05i0.13 1 


















BPL 306 




19. 09 9 


16.16±0.04 


15.53±0.02 


15.15±0.02 9 


y 


y 


y 










BFL 316 




18. 30 9 


15.97±0.04 


15.41±0.02 




y 


y 


y 




















14.93±0.04 
















BPL 327 


IPMBD 11 


17.94 9 


15.52±0.03 


14.99±0.02 


14.60±0.02^ 


y 


y 


y 




y 


y g 




PPL 1 


Roque 15 


18.21* 


15.37±0.08* 


14.82±0.06* 


14.34i0.05 1 


y 


y 


y 


b 


y 




y 


PPL 15 


IPMBD 23 


18.24' 


15.34* 


14.77* 


14.4P 


y 


y 


y 


b 


y 


y 


y 




NPL 35 
























CFHT-PL-23 




19.33* 


16.38±0.08* 


15.79±0.08* 


15.25i0.05 1 


y 


y 


y 


b 








Roque 5 




20.19* 


16.89±0.14* 


16.33±0.15* 


15.74i0.10 1 


y 


y 


y 










Roque 8 




19.55' 


16.68±0.11* 


16.14±0.10* 


15.57zb0.06 1 


y 


y 


y 










Roque 25 




21.80* 


17.7±0.3* 


17.0±0.3* 


16.31zb0.19 1 


y 


y 


y 




y 






Roque 30 




20.82* 


17.49±0.12* 


16.77i0.ll* 


16.14ib0.07 1 


y 


y 


y 










Roque 33 


NPL 40 


20.44* 


17.07±0.12* 


16.62±0.18* 


16.02i0.08 1 


y 


y 


? 




y 


y j 


y 


Roque 36 




20.25* 


17.17±0.20* 


16.58±0.12* 


16.09zb0.08 1 


y 


y 


y 
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Table A2. Photometry (c2-5) and membership (c6-8) of Pleiades low-mass star candidates (higher mass BPL candidates). Membership 
is from the IK CMD, the JK CMD and the JHK diagram. J- H- and K-band photometry is from 2MASS, and has been transformed into 
the MKO system. 



BPL 


Ic 


Jm KO 




Kjv/if 


IK 


JK 


JHK 


3 


14.81 


i q 9^+0 03 


1 2 74+0 03 

li. l1_l_U.UO 


1 9 3fi+n ri3 


y 


y 


y 


4 


14.17 


i o fi7-U0 03 


1 9 1 94-0 09 

_1_Z . l^XU. uz 


1 1 774-0 03 

1 1 . 1 ( XU. UO 


y 


y 


y 


5 


14.25 


1 9 fifi+0 03 


1 2 1 7+n 

li.l 1 _I_U. uo 


1 1 so+0 03 


y 


y 


y 




15.74 


1 4 ni +n 03 

J — l_U . UO 


1J . oo_i_u . uo 


1 q nQ+n ri4 


y 


y 


y 


7 


16.28 


1 A 44+0 03 


1 q 094-O 03 
10 . yzxu. uo 


1 q A^+0 04 

10.10IU.U4: 


y 


y 


y 


16 


17.02 


1 5 20+0 04 

lij . iiU_l_U. U1 


1 a (S74-0 O^i 


1 4 94+0 Ofi 


N 


y 


y 


18 


14.39 


1 O QQ-UO 03 

iz . oyiu. uo 


1 9 3^4-0 03 

IZ . OOXU.UO 


1 1 074-O 03 
ii.y 1 xu. uo 


y 


y 


y 


19 


16.39 


1 A fiO+0 03 


1 a 0Q4-0 OA 

_i_1. Uy XI U. Ul 


1 q 7^4-0 04 
10. 1 uxu.ui 


y 


y 


y 


90 

ZU 


14 81 


1 Q 91 4-0 03 


1 9 704-0 03 

IZ. f UXIU.UO 


1 9 9«4-0 03 

IZ.ZOxU.UO 


y 


y 


y 


26 


17.00 


i a fi74-0 03 

11. Ul IU.UO 


1 a 1 R4-0 OA 

XL 1UXU. Ul 


1 q 7^4-0 O^t 

10. 1 UXU.UO 


y 


y 


y 


30 


15.90 


i a 1 c;4-0 03 


1 q c;q4-o 03 

10 . OuXU. uo 


1 q 994-O 04 

J.O. ZZiXZU. U4: 


y 


y 


y 


34 


13.72 


i 9 m 4-0 09 

l^.Ul XU. UZ 


1 1 Ac;4-0 03 

1 1 . liJXU.UO 


1 1 1 0+0 03 

1 1 . 1UXU. uo 


y 


y 


y 


38 


16.12 


i 4 974-0 03 

11. Zf XU.UO 


1 q 7Q4-O OA 

IO. / OXU.U4 


1 q A9+0 04 

10.4:ZrXU. U4: 


y 


y 


y 


43 


15.81 


1411 J_n riQ 

14:. 1 1 XU. UO 


1 q c;(S4-0 03 

IO . JUXU.UO 


13 1 8+0 04 

IO. IOIU.UI 


y 


y 


y 


44 


16.09 


1 a 304-0 03 


1 q 7R4-0 03 
10. / uxu.uo 


1 q 374-O 04 

IO.O t XU.U4: 


y 


y 


y 


^fi 


1 4 on 
14:. yu 


1 q 9Q4-0 03 
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1 9 7A4-0 09 
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IZ.'iUxU.UO 


y 


y 


y 


58 


17.69 
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y 


y 


y 


ui 


14:. OU 


1 Q 1 1 4-0 09 

IO. ilXU.UZ 


1 9 c:o_Ln 03 


1 9 934-0 03 
IZ.ZOxU.UO 


y 


y 


y 


68 


14.33 


1 9 Qccj-n 03 
1Z . ooxu . uo 


1 9 3A4-0 03 

li. 04:XIU. UO 


1 9 03+0 04 
1z.uoxu.u1 


y 


y 


y 
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uy 
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lU.OlxU. 10 
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1U. lUXIU.ZU 
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N 


y 


y 


70 


15.40 


1 q 334-O 09 
IO . ooxu. uz 


1 9 7(S4-0 03 


1 9 a^+0 03 
1z.101u.uo 


y 


y 


y 


71 


15 14 


1 q 1 Q4-n 09 

IO. lOXU.UZ 


1 9 c:o_Ln 03 


1 9 934-O 03 

IZ.ZOxU.UO 


y 


y 


y 


72 


16.18 


1 A OA+0 03 


1 q c;94-0 03 
10 . uuxu.uo 


1 q 90+O 03 

lO.ZUIU.UO 


y 


y 


y 


73 
f O 


1 ^ OQ 


1 Q OS4-0 09 
io.uoznu.uz 


1 9 A04-0 03 
iz.uuznu.uo 


1991 4-0 03 

1 Z . ZlxU.UO 


y 


y 


y 


74 


17.44 


1 cr 1 Q4-0 O 1 ^ 

10 . lyiu.uo 


1 a fiS4-0 Ofi 

l'i.UOXU.UU 


1 a 934-0 Ofi 

11.ZOIU.UU 


y 


y 


y 


7^ 


1 fi 98 


1 A 034-0 03 

14:.UOxU.UO 


1 3 Afi4-0 03 
io.4:Ux;u.uo 


1 q 0fi+0 03 

lO.UUXU.UO 


y 


y 


y 


84 


17.21 


1 a Qi 4-0 OA 

11.01 XU. Ul 


1 a 9Q4-0 OA 

I'i.ZiCJXU. U4: 


1 q ec;4-o 04 

lO.OOIU.UI 
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y 


y 


oy 


io.oy 
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IZ.UlxU.UO 


1 1 Afi4-0 03 
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y 


y 


y 


91 


14.74 


1 9 734-0 03 
IZ. f oxu.uo 


1 9 01 4-o 03 

11,^.1 XI u. uo 


1 1 QQ4-0 03 
1 1 . ooxu. uo 


y 


y 


y 


92 


15.44 


1 q 394-O 03 
10 . O^XIU. uo 


1 9 Q34-O 03 

li. OOXU.UO 


1 9 A4+0 03 
1z.11xu.uo 


y 


y 


y 


93 


14.63 


1 9 71 4-n nQ 


1 9 1 24-0 03 


1 1 01 4-0 03 


y 


y 


y 


QA 


1 ^ 9fi 


1 Q c;fi4-0 03 
IO. JoxU.UO 


1 q 0S4-0 03 
10.uon1u.uo 


1 9 74+0 03 

1Z. / IxU.UO 


y 


y 


y 


95 


14.02 


1 9 e:i 4-n no 

IZ . O 1XU.UO 


1 1 Q(S4-o 03 
1 1 . yuxu.uo 


1 1 (St;4-o 03 

1 1 . UOIU.UO 


y 


y 


y 


Qfi 


1 fi 38 


1 A 0^4-0 03 

14:.UUXU.UO 


1 Q c;i 4-0 OA 

IO. 01XU.U4: 


1 q 1 Q4-0 03 

IO. IOIU.UO 


y 


y 


y 


97 


14.55 


1 9 Q94-0 03 
IZ . ozxu. uo 


1 9 394-0 03 


1 1 Qfi4-0 03 
1 1 . yuxu.uo 


y 


y 


y 




1 ^ 89 
lO.OZ 


1 q Qp;4-0 03 
10. yuxu.uo 


1 3Q4-0 03 
10.oyxu.uo 


13 11 +0 04 
10. iixu.ui 


y 


y 


y 


99 


15.30 


1 q fi0+0 03 
10 . uuxu.uo 


1 q 0S4-0 03 
10 . UO X u. uo 


1 9 744-0 03 

1Z. 1 1XU.UO 


y 


y 


y 


1 01 

lul 


1 fi 70 
1U. I u 


1 A 7A+0 OA 

14:. ( 4:XIU.U4: 


1 A 1 Q4-0 0^ 
14:. lyxu.uo 


1 q Q94-0 Ofi 

lO.OZxU.UU 


y 


y 


y 


102 


15 11 


1 3 40+0 03 
10.4:UXIU.U0 


1 9 Q74-0 03 

1Z.O 1 IU.UO 


1 9 c;c;4-0 03 

1Z.OOXU.UO 


y 


y 


y 


103 


14.57 


1 9 Q04-0 03 
iz . yuxu.uo 


1 9 384-0 03 
1^. ooxu . uo 


1 9 o^+O 03 


y 


y 


y 


104 


16.26 


1 a 0Q4-0 03 


1 q c;a4-0 OA 

IO . 04:XU. U4: 


13 1 Q+0 03 

lO.l U^.\J'UO 


y 


y 


y 


1 0^ 


1 ^ AQ 


1 q 334-0 03 

lO.OOxU.UO 


1 9 8A4-0 03 

1Z.04:XU.UO 


1 9 A3+0 03 

1Z.1iOxU.UO 


y 


y 


y 


106 


17.44 


1 k 994-0 Ofi 

IO .ZZXU. UU 


1 a fi94-0 07 

14:. UZXU.U 1 


141 c:4-0 07 

11. IOXU. u t 


y 


y 


y 


107 


16.18 


1 A 04+0 03 


1 q c;a4-0 OA 

IO . l)1IU. U4: 


13 1 fi+0 04 
10. luxu.ui 


y 


y 


y 


109 


16.02 


1 a 1 C4-0 03 

J_1. lUXU. UO 


1 q 70+O OA 


1 3 9S+0 O'S 

J-O. iiO_l_U. uo 


y 


y 


y 


110 


16.70 


14.40±0.03 


13.83±0.04 


13.47±0.05 


y 


y 


y 


111 


15.98 


14.19±0.03 


13.62±0.04 


13.23±0.04 


y 


y 


y 


112 


17.30 


14.98±0.05 


14.37±0.06 


14.21±0.08 


y 


y 


y 


113 


13.73 


12.02±0.02 


11.45±0.03 


11.12±0.02 


y 


y 


y 


114 


13.85 


12.34±0.02 


11.82±0.03 


11.48±0.03 


y 


y 


y 


115 


17.00 


14.90±0.05 


14.37±0.05 


13.98±0.06 


y 


y 


y 


116 


15.79 


13.75±0.03 


13.27±0.03 


12.87±0.03 


y 


y 


y 


117 


14.64 


13.00±0.03 


12.49±0.03 


12.12±0.03 


y 


y 


y 


118 


14.82 


13.07±0.02 


12.53±0.03 


12.17±0.03 


y 


y 


y 


119 


15.10 


13.23±0.03 


12.67±0.03 


12.29±0.04 


y 


y 


y 


120 


15.06 


13.11±0.03 


12.59±0.03 


12.23±0.03 


y 


y 


y 


121 


15.03 


13.20±0.03 


12.67±0.03 


12.33±0.03 


y 


y 


y 
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BPL 


lc 


Jmko 


H_MJfO 


Kjv/if 


IK 


JK 


JHK 


122 


17.13 


14.96±0.05 


14.46+0.06 


14.04+0.07 


y 


y 


v 
J 


123 


17.32 


1 5 74+0 08 

i '-J . I I kj . w yj 


14 98+0 09 


1 4 90+0 1 5 


N 


y 


N 


125 


15.80 


1 3 98+0 03 

_L U . C/KJ _1_ U . W U 


1 3 40+0 04 


1 3 08+0 03 


y 


y 


y 


126 


14.34 


12 74+0 02 


12 18+0 03 


1 1 87+0 03 

J- -L • C 1 V7 * UU 


y 


y 


y 


127 


17.71 


1 5 37+0 06 

UiO 1 KJ . yjyj 


14 84+0 07 


14 32+0 08 


y 


y 


y 


128 


16.45 


1 4 44+n 03 


1 q qi +0 04 


1 q R4+0 05 


y 


y 


y 


129 


16.23 


1 4 23+0 03 


1 q 7R+0 04 


1 q 4q+o 04 


y 


y 


y 


130 


17.27 


i 4 ai +n (14 


1 4 27+0 05 


1 q 07+0 06 


y 


y 


y 


131 


16.12 


14 00+0 03 


1 3 55+0 04 


1 3 1 8+0 04 

J. U • J- yj 1 v7 * utt 


y 


y 


y 


133 


14.02 


12.34+0.03 


11.80+0.03 


11.45+0.03 






y 


134 


15.60 


1 3 62+0 02 


13 12+0 03 


1 2 75+0 03 

J- — ■ 1 \J 1 w * yj * j 


y 


y 


y 


135 


15.22 


1 3 21 +0 04 


12 73+0 03 

-L ■ J O yj . yj >.J 


1 2 36+0 03 

1 • ' 7 y j 1 y ' - v ' ■ ) 


y 


y 


y 


136 


15.29 


1 3 37+0 03 

_i_u.ui kj . w u 


12 84+0 03 


1 2 46+0 04 


y 


y 


y 


138 


15.09 


1 3 14+0 02 

_i_ u . -L I yj . w ±j 


12 65+0 03 


1 2 29+0 03 

J_ — ■ £j Zf 1 KJ * KJ *-J 


y 


y 


y 


139 


16.23 


14 07+0 03 


1 3 53+0 03 


1 3 1 7+0 04 

J_U» J- 1 1 U-U^ 


y 


y 


y 


140 


16.44 


1 4 22+0 04 


1 q 71 +n 04 


1 q q.2+0 04 


y 


y 


y 


141 


14.27 


i 2 fi4+n 03 

. Urt_l_U. UU 


1 2 1 q+o 03 


1 1 7S+0 03 


y 


y 


y 


143 


15.70 


i q S4+n 03 

J_U . Ort_l_U . UiJ 


1 q q.4+0 03 

iij . Ort_l_L/. KJiJ 


1 q 00+0 03 

J-tJ. > ' 1 ' ' . KJ' J 


y 


y 


y 


144 


14.06 


i 2 q.0+0 02 

li . 'J KJ _L KJ . Ui 


1 1 7^+0 03 

11.1 UJ.U . wo 


1 1 41 _|_n 03 


y 


y 


y 


145 


15.34 


1 3 68+0 03 


13 17+0 03 

-L y J . -L 1 _l_ U . W >.J 


1 2 86+0 04 

1 ^- - 1 1 yj 1 v / . vji 


y 


y 


y 


146 


16.77 


14 97+0 05 

_Lrr . C 1 _1_ U . UU 


14 40+0 06 


14 21+0 07 


N 


y 


y 


147 


15.50 


1 3 69+0 03 

J-UtUC _1_ U . UU 


1 3 22+0 03 


1 2 87+0 03 

j_ — ■ <■ ? 1 1 v 7 1 yj*_j 


y 


y 


y 


148 


14.93 


i 2 Q7+n 04 

li.J 1 _l_ U . Urt 


1 2 42+0 03 


1 2 07+0 03 


y 


y 


y 


149 


15.45 


i q R7+n 03 

-LU.Ul -1— KJ . KJ i J 


1 q 1 n+n nq, 

-LO . -LU_I_LJ. KJyJ 


1 2 7^+0 04 


y 


y 


y 


150 


15.46 


13.64+0.03 


13.12+0.03 


12.74+0.03 


v 

y 


v 

y 


v 

J 


151 


14.04 


1 2 55+0 03 

i _ . u ' ) v> .y > • > 


12 03+0 03 


1 1 68+0 03 

J_ _L m KJ O 1 yj * UU 


y 


y 


y 


153 


14.76 


i q 03+0 03 

_L U . U U _J_ U . U U 


1 2 4Q+0 03 


12 1 q+o 03 


y 


y 


y 


154 


16.94 


1 4 72+0 04 


1 4 20+0 05 

it.iUiu . yj'j 


1 q 7S+0 05 

J-U. 1 ' ' 1 ' ' . KJyJ 


y 


y 


y 


155 


15.52 


1 3 49+0 03 

_l u . rt c _i_ u . uu 


12 88+0 03 


1 2 54+0 03 

1 ^- - ' ) 1 1 v / . uu 


y 


y 


y 


156 


15.30 


13 41+0 03 

i- U i T: J — 1_ U . UU 


12 91+0 03 

j. _ . ■ 1 \ — 1_ yj . yj >.j 


1 2 56+0 04 

j_ ^ . \7 yj 1 v ) . utt 


y 


y 


y 


157 


14.78 


i 2 q^+0 03 


1 2 42+0 03 


1 2 07+0 03 
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y 


158 


15.27 


i q 4«+o 03 


1 q ni +0 03 

-L i J . yj J — 1_ KJ . KJ i J 


1 2 R2+0 03 

J- . UZ/-LU. UU 
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y 


159 


13.78 


i 2 q.0+0 03 

li . JU.LU . UU 


1 1 77+0 03 

-L -L . ( i —1— KJ . KJ i J 


1 1 4^+0 03 

J. -L .rtU—LU. UU 
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y 


y 


160 


15.06 


12 96+0 03 

_1_ _ . CU _1_ U . UU 


12 47+0 04 


12 12+0 03 

J_ — ■ -L 1 yj * UU 


y 


y 


y 


161 


15.78 


13.72+0.04 


13.20+0.04 


12.87+0.03 


v 
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v 
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v 


162 


16.07 


14 15+0 04 

iTtt lu _1_ U . UT: 


1 3 62+0 05 
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1 3 25+0 05 

1 -J - 1 v, 7 . uu 
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y 
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16.38 


1 4 43+0 04 


1 q qfi+o 05 
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y 


165 


14.04 


1 1 81 +0 03 

iiiUJ — 1_ u . uu 


1 1 37+0 04 


10 97+0 03 
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y 
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y 


166 


14.09 


1 2 35+0 04 


1 1 82+0 03 


1 1 48+0 03 
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y 


y 


167 


13.74 
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1 1 55+0 05 

J- J. . yj ' J \ KJ . KJ ' J 
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y 


y 


y 


169 


17.60 


15 32+0 06 

i ' ) . - ) — u . uu 


14 71+0 07 


14 43+0 08 


y 


y 


y 


170 


15.68 


1 Q QQJ_n HA 


10.0 1 xu. uu 


1 q7-un 
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y 


y 


171 


15.43 
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y 
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14.16 
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13.18±0.04 
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y 


y 


y 


184 


14.21 


12.68+0.03 
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y 
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y 
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y 


y 


y 


285 


14.97 


13.36+0.03 
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14.46+0.06 
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Table A3. Photometry (c2-5), membership (c6-8) and binarity (c9) of RlZ-Pr Praesepe candidates. Membership is from the IK CMD, 
the JK CMD and the JHK diagram. Refs; (a) Pinfield ct al. (1997), (b) Hodgkin et al. (1999). The t is also IZ-Pr 40. The * is also IZ-Pr 
86. n-d = non-detection 



RlZ-Pr \ a c 3 MKO Kmko K mko IK JK JHK Bin 



1 


20.14 


17.06i0.05 


16.47i0.04 


16.14i0.08 6 


y 
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y 


2t 
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15.67±0.05 


15.11±0.03 
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y 


y 
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3 
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y 


N 


N 


4 
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15.56±0.03 
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y 


y 


y 


5 
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N 


N 


6 
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21.24 


- 


- 


18.56±0.29 
19.15i0.05 6 
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Table A4. Photometry (c2-5), membership (c6-8) and binarity (c9) of IZ-Pr Praesepe candidates. Membership is from the IK CMD, 
the JK CMD and the JHK diagram. Refs; (a) Pinfield etal (2000), (b) 2MASS. The t is also FJZ-Pr 2. The * is also RIZ-PR 21. 
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